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Gluino pair production in an EGRET motivated super symmetry scenario has been 
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studied. The analysis shows, that this special scenario could be proven with the first 
LHC data, and assuming a higher integrated luminosity the analysis would be able 
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to cover a large region of the minimal supergravity (mSUGRA) parameter space. 
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I. INTRODUCTION 

The Minimal Super symmetric Standard Model (MSSM) is an extension of the standard 
model including a symmetry between fermions and bosons. Such models introduce fermion 
partners to all standard model particles and bosonic ones to all fermions (SUSY particles). 
Assuming, that typical SUSY masses are located at the TeV scale would make the unifi- 
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cation of the coupling constants easily accessible Additionally, in R-parity conserving 
models the neutralino would provide a candidate for dark matter, a so-called WIMP (weakly 
interacting massive particle) [2]. 
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Due to the fact, that supersymmetry (SUSY) has not been discovered yet, high masses 
for SUSY particles can be achieved by applying Soft Symmetry breaking mechanisms 
These symmetry breaking mechanisms results in general in more than 100 free parameters, 
which are reduced to five in minimal supergravity (mSUGRA) models: the common scalar 
(fermion) mass parameters mo (miw) at the GUT scale, the sign of the Higgs mass param- 
eter sign(//), the trilinear coupling A Q , and the ratio between the two neutral Higgs vacuum 
expectation values of the two Higgs doubletstan (3 [31]. Thus the SUSY parameter set is given 
by 

{p} = {sign(/i),TOo,mi/ 2 , A ,tan/?}. (1) 
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Due to the fact, that the masses of fermions, hence their cross sections and decay channels, 
strongly depend on the chosen mSUGRA parameters an analysis cannot cover the whole 
mSUGRA parameter space. Thus a reasonable region for the analysis should be chosen. 
This choice can be motivated by the EGRET experiment, with an excess of high energetic 
galactic gamma radiation has been observed. This excess can be explained by dark matter 
annihilation with neutralino mass in the range 50 to 70 GeV j^. This condition corresponds 
to a special regions in the mSUGRA parameter space: The region at low fermion masses, 
which corresponds to low m^j. Here the mass of the neutralino is given by pa 0.4mi/2, 
while the mass of gluinos is rrig pa 2.7mi/2- In this region the masses of x® and Xi are about 
one fourth of the gluino mass. The value for the dark matter density given by these models 
would be consistent with the results of WMAP [2J. The regions, which are preferred by 
experimental and theoretical results are shown in Fig. [1] [4|. 

Another region, which is consistent with the results of WMAP is the focus point region 

n 

(FPR) |5J. This region is a narrow band along the region, which is forbidden by electroweak 
symmetry breaking. In this region the masses of the two lightest neutralinos get very close 
to each other, which makes them independent of m^. The mass of the gluino is expected 
to be about one order of magnitude larger than the ones ofj^ and Xi- 
Thus the cross sections for gaugino production are expected to be higher compared with 
the cross sections of sfermion production. The production channels with the largest cross 
sections in this scenario are expected to be gluino pair gg production and direct production 
of neutralino chargino X2X1 pairs, which can be detected as a trilepton signature jf|. 
The relation between the two production channels can be determined from the mass 
hierarchies, as they have been discussed before. In the low m^j region the cross section for 
gluino pair production is expected to be one order of magnitude higher than the one for 
neutralino chargino pairs, which is shown in Fig. [2j while in the FPR neutralio-chargino 
pair production is preferred. 

The trilepton signature is hard to separate from the background. Thus it will be harder 
to detect than gluino pair production. Additional, the decay cascades of gluino production 
are expected to be longer, thus more parameters can be determined from the same data. 
Gluino pair production is in both regions the gluino production channel with the highest 
cross section. The Feynman diagrams, which contribute to gluino pair production are given 
in Fig. [3] |3(. In the low mi/2 region the production cross sections are large compared to 
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other regions (in the 10 pb range) caused by the low mass of the gluino. In the focus point 
region the gluinos are much heavier, thus the production cross sections vary between pb 
and fb. 

In both regions the gluino is lighter than any squark. This results in three-body decays 
via virtual squarks. The products continue decaying into lighter SUSY particles until the 
LSP is reached. Thus the length of the decay chain depends on the mass of the gluinos and 
the gauginos. The branching ratios of the decays to bottom and top quarks are strongly 
related to the initial values of the SUSY breaking. For the low fermion mass region the 
decay to bottom quarks is strongly preferred. Thus the main decay channel is the decay 
into a neutralino and a bottom pair and the decay into a bottom quark. Caused by the low 
mass of the gluino the decay cascades appear not to be very long. In the lower part of the 
focus point region also the decay to top quarks contributes and in the higher part also to 
channel including a bottom, a top and a chargino is open. The main contributing Feynman 
diagrams for gluino decays are plotted in Fig. HI Two examples for cascades at LM9 and in 
the FPR are given in Fig. [5j 



II. ANALYSIS 



The analysis has been done at Monte Carlo level, i.e. excluding detector simulation and 



reconstruction, 
table is given in 



?he hardware based LI trigger from CMS has been included. Its trigger 
7| and the resulting trigger efficiency of the signal and background events 
are all mentioned in table [H 

The event topology for gluino pair production, which can be assumed from the possible 
types of cascade decays, is given by: 

Ef ss high 

™jets > 4 

M eS = E™ s + J2PT,high (2) 

jets 

^-leptons • 

The missing transverse energy (E™ 1SS ) is expected to be high, because the neutralino is 
expected to escape from the detector. The third condition is caused by the high sfermion 
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masses. In the gluino decay cascades leptons can only be created in decays of neutralinos 



or charginos. In order to conserve colour charge minimal four 



ets have to be produced per 



event. Additionally the effective mass is expected to be high js|. That is caused by the fact 
that SUSY particles are heavy, thus the MET and the transverse momenta of all jets in the 
events are expected to be high. Cased by the fact, that the gluino decays into electroweak 
gauginos, also leptons can appear in the event topology. 

Possible standard model backgrounds for the production of heavy SUSY particles are heavy 
quark pair production (tt, bb and cc) and heavy vector boson production (W and Z plus jets). 
The heavy quark pair production contributes to the missing transverse energy and to the 
number of jets. Thus the worst background is expected to be tt, because it is dealing with 
the heaviest particles, but its cross section is low (less than 1 nb) compared with the other 
backgrounds. Many orders of magnitude larger cross sections are expected from vector boson 
production (in the 50 nb range). That is, what makes these backgrounds very dangerous. 
They contribute to the missing transverse energy, to the number of jets and from the decays 
of the vector bosons also to the number of leptons. 

The mSUGRA parameter point our analysis has been tuned to the CMS benchmark point 
LM9. It is located in the low fermion mass region and its parameter set is given by 

{p} = {+1, 1450 GeV, 175 GeV, 0, 50} (3) 

The cross section at this parameter point can be calculated using Prospino [9], and it is 
20.8 pb in leading order and 36.7 pb in the next to leading order approximation. 
The given event topologies for the low fermion mass point LM9 (definition given in ([3])) are 
plotted in Fig. [6] and [3 In reality the topologies of the events are distorted by jet and 
MET misreconstruction. There it can be seen, that the MET is an important variable to 
distinguish between signal and background events. 

The separation between signal and background was optimised using neural networks (Neu- 
robayes [lQ|). To check the consistency of the neural network the analysis has also been 
done as cut analysis. 

In a neural network the variables get different significances. The more the variables con- 
tribute to the output, the higher their significances. Variables, which have been used to train 
the neural networks and their significances are given in table HI1 The two most significant 
variables are M e g and E™ 1SS . Therefore one can try to separate signal from background in a 
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2D-scattering plot of M e g and E™ 1SS as shown in Fig. [HJ The event numbers are shown in 
the second row of table IHIl 

A neural network is sensitive to statistics, thus it cannot be trained to a sample, where 
the number of background events is much higher than the one of signal events. Because of 
the high background cross section precuts have to be applied to the data. The following 
preselection cuts have been done: 

M eff > 500 GeV 

(4) 

E™ ss > 150 GeV. 

The resulting event numbers are written in the third row of table IHIl After applying 
the preselection the numbers of background events are in the same order of magnitude as 
the number of signal events, which provides perfect conditions for the training of a neural 
network. 

To avoid a loss of information three different neural networks have been trained to different 
groups of backgrounds. Here backgrounds with high similarities have been taken as one 
group: WJets and ZJets and bb and cc. ti has been chosen as separate group. 
After the training the problem can be reduced to a problem in a three dimensional variable 
space of the network outputs. Because the neural networks have been trained separately, 
the outputs have to be normalised to the different cross sections of the backgrounds. 
To do the final event selection on the strongly correlated outputs of the three neural networks 
another neural network was trained to the outputs. Thus the problem is in the end reduced 
to a one-dimensional problem. An overview of the neural network analysis is given in Fig. 
[9j The obtained one dimensional problem can easily be solved by significance optimisation. 
The resulting event numbers are listed in table IHIl 

To proof a discovery the significance, which is for a Gaussian distribution defined as 

S 



a 



(5) 



where S is the number of signal events and Bi is the number of background events for the 
background i, has to be more than five. 

With the given results the significance at LM9 can be calculated to be a m 60. From 
comparison of the results for the 2D cut method and the neural networks it can be seen 
that the number of background events is comparable for both cases. What makes the neural 
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network analysis better is the signal efficiency. In the neural network analysis about double 
the number of events is selected correctly. The fact, that the numbers of background events 
are in the same order of magnitude after both selections, validates the consistency of the 
neural network. 

The analysis has been tuned to one special parameter point. Now it should be checked, how 
general the analysis is and if it could be applied to other points in mSUGRA parameter 
space. Thus a scan in the m Q -mi/2 plane has been done at fixed tan (3=50. The results are 
plotted in Fig. [TJJJ 

The analysis can be applied for parameter points surrounding LM9. It is sensitive to the 
gaugino mass parameter mm and less sensitive to the scalar mass parameter mo, which is 
not surprising, because the gluino is a spin 1/2 particle. As shown in Fig. \TU\ the fermion 
mass parameters up to 500 GeV and scalar mass parameters up to 3000 GeV are within the 
reach for an integrated luminosity of 30 fb _1 , which corresponds to three years LHC running 
time. 

In the focus point region (description given in |5j) a more detailed parameter scan has been 
done. The focus point region is interesting for the comparison, due to the fact that in this 
region the event topology differs from the one at LM9. Here the decay to bottom quarks 
is not strongly preferred since the scalar mass parameter becomes high enough to give the 
possibility for a proper coupling to the top quark. The results are shown in Fig. [TTJ 
The analysis works for some lower points in the focus point region, but it does not in the 
upper area. Here the event topologies get too different from the ones at LM9 

III. CONCLUSION 

The EGRET experiment motivates the CMS scenario LM9 ([3]) and the focus point region. 
The Monte Carlo analysis shows that an mSUGRA the significance for a discovery of an LM9 
scenario (tan/? = 50) would be about 60 with an integrated luminosity of 1 fb" 1 , which is 
reached after one year LHC runtime. With an integrated luminosity of 30 fb _1 the analysis 
would cover the low fermion mass region of (m , mi/2) = (200 — 2800, 100 — 400). 
The next steps, are including the detector simulation, which is expected to decrease the 
significance. Additionally, the event selection should be optimised with real data, as soon 
as available. 



[1] W. de Boer and C. Sander, Physics Letters B 585, 276 (2003). 

[2] W. de Boer, M. Herold, C. Sander, V.Zhukov, A. Gladyshev, and D. Kazakov (2007), 

arXiv: astro-ph /0408272 . 

[3] S. P. Martin (2006), arXiv:9709356 [hep-ph]. 

[4] W. de Boer, CERN Cour. 45, 17 (2005). 

[5] H. Baer, T. Krupovniackas, S. Profumo, and P. Ullino, JHEP 0510, 020 (2005). 

[6] W. de Boer, I. Gebauer, M. Niegel, C. Sander, M. Weber, V. Zhukov, and K. Mazumdar 

(2006), cERN-CMS-NOTE-2006-113. 

[7] C. Collaboration, J. Phys. G: Nucl. Part. Phys. 34, 995 (2006). 

[8] M. Chiorboli (2007), arXiv:0710.4787 [hep-ph]. 

[9] W. Beenakker, R. Hopker, M. Spira, and P. M. Zerwas, Nucl. Phys. B 492, 51 (1997). 

[10] M. Feindt (2004), arXiv:physics/0402093vl. 

[11] T. Sjostrand, L. Lonnblad, and S. Mrenna (2001), arXiv:hep-ph/0108264. 



8 



Figure 1: Constraints on the mSUGRA parameter space for tan/3 = 50. The region, which is 
allowed by the results of WMAP are assigned in blue. The EGRET data require low m-1/2 values, 
which implies values of m-o in the tera-electron-volt range, if all constraints from the Higgs limit 
(rrih), g — 2 and the b — > S7 branching ratio are considered (limits exclude regions to the left). The 
CMS point LM9 (mo = 1450, tti 1 / 2 = 175) is located in the region which is preferred by theory 



and experiments 
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Figure 2: Cross sections of SUSY processes at the CMS parameter point LM9 normalised to the 
sum of all SUSY cross sections versus the centre of mass energy of the collider. At the Tevatron 
(about 1000 GeV) neutralino-chargino pair production x^X° would be the dominant process, while 
the dominant process at the LHC (14000 GeV) will be gluino pair production (gg). All processes, 
which are not plotted (like squark pair production qq), contribute less than 10 to the SUSY cross 
section. 
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Figure 3: Feynman diagrams for gluino pair production. Gluino pairs can be produced from a 
gluon directly (like plotted in the first line) or via exchange of a gluino or a squark (like plotted 
in the second and third line). In the LHC production from quark pairs (right row) as well as 
production from gluon pairs (left row) contribute to the production cross section. 
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Figure 4: Feynman diagrams for gluino decays at LM9 and in the focus point region. For LM9 the 
first diagram is strongly preferred compared to the other ones. In the focus point region, at low 
scalar masses the branching ratio of the second digram gets two times larger than the one of the 
first diagram. At higher fermion masses (thus also higher scalar masses) the channels, which are 
plotted in the lower line, are open. 
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Figure 5: Two examples for decay cascades at LM9 and in the FPR. At LM9 the decays to top 
quarks are strongly suppressed, because the scalar mass parameter is very low. The neutralino re- 
sults in missing transverse energy, because it is electrically neutral and weakly interacting. Leptons 
contribute to the lepton number and jets to the jetmultiplicity as well as during their transverse 
momenta to the effective mass. 
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Figure 6: Number of jets versus number of leptons. The event numbers have been normalised 
to one. Due to jet missreconstruction the number of jets gets similar to the one of gluino pair 
production. As result the topology of ti production is very similar to gluino pair production. 
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Figure 7: MET versus number of jets. The event numbers have been normalised to one. The MET is 
very different for ti events and gluino pair production. That is because the lightest supersymmetric 
particle, which is very heavy, is escaping from the detector. The MET can be used very well to 
distinguish between signal and background. 
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Figure 8: Scattering plots including two dimensional cut function. The scattering plots are nor- 
malised to one, thus the ones of the backgrounds include a lot more events. The black lines assign 
the two dimensional cuts, which have been applied to check the consistency of the neural network 
analysis. The grey lines show the precuts, which have been done to prepare the data samples for 
the neural network training. 
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Figure 9: Structure of the neural network analysis. The distributions of the backgrounds are not all 
similar to each other. To avoid a loss of information the analysis has been divided into two parts: 
The training to different background samples and the training to the outputs of the different neural 
networks. Twelve variables are given as input for the first step neural networks: I?™ 1SS , nLeptons; 
"Jetsj ^B-Jets, M e g, cos (C)> an d Pt and rj of the three leading jets. The final neural network is 
trained to the outputs of the of the previous three neural networks. As last step a one dimensional 
cut is done on the final neural network output, which looks well separated after applying the final 
neural network. 
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Figure 10: Significance scan in the mSUGRA mo-m^j plane at tan/3 = 50. for an integrated 
luminosity of 30 fb _1 . The discovery reach is up to = 500 GeV. The lower right corner is 

exculded by EWSB and the upper left part is not in accord with WMAP data. 
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Figure 11: As in Fig. [TUl but for a finer scan in mo and m^ 2 . Here just for the lower points the 
significance is high enough to make a discovery. 
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Table I: Evolution of event numbers during the analysis for 1 fb 1 . First the event numbers without 
any analysis tool are plotted. The px (definition given in l_l|]) cuts chosen for the generation are 



the following: 50 GeV for ZJets, 85 GeV for WJets, and 180 GeV for bb and cc. In the first step 
the CMS LI trigger is included, where WJets and ZJets already get reduced very well, while 87% 
of the signal events pass the trigger. Backgrounds that also pass the trigger nearly untouched are 
tt and bb. After applying the preselection cuts the number of background events is in the same 
order of magnitude like the number of signal events. Also the number of tt is just about four times 
bigger than the number of signal events. After applying the final selection the final significance 
can be calculated to be a = 54. 
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Table II: Significances of the variables used for the neural network training. The three jets chosen 
are the three hardest jets in the event and ( is the angle between the missing transverse momentum 
and the transverse momentum of the hardest jet. 
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Table III: Results of the two dimensional cuts and the neural network. The outputs of the neural 
network analysis and the 2D analysis are in the same order of magnitude, which validates the 
neural network analysis. For the neural network analysis the signal efficiency is higher and the ti 
and the bb as well as the ZJets selection is better, while the cc and WJets selections are a bit worse 
in the neural network analysis. That is caused by the fact, that these backgrounds are trained 
together with other backgrounds. So if one background separation gets better the other will get a 
bit worse. But in case of our analysis the avail is larger than the loss. 
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